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The Space Interferometry Mission (SIM) is a space-based long-baseline optical interferometer for precision astrom- 
etry. One of the primary objectives of the SIM instrument is to accurately determine the directions to a grid of 



stars, together with their proper motions and parallaxes, improving a priori knowledge by nearly three orders of 

I , 

magnitude. The basic astrometric observable of the instrument is the pathlength delay, a measurement made by a 

Ctf ■ 

Oh' combination of internal metrology measurements that determine the distance the starlight travels through the two 
arms of the interferometer and a measurement of the white light stellar fringe to find the point of equal pathlength. 

O i' the interferometer baseline vector is not stationary over this time period, as its absolute length and orientation 



Because this operation requires a non-negligible integration time to accurately measure the stellar fringe position, 



are time-varying. This conflicts with the consistency condition necessary for extracting the astrometric parameters 

>'• 

which requires a stationary baseline vector. This paper addresses how the time-varying baseline is "regularized" so 

o 



that it may act as a single baseline vector for multiple stars, and thereby establishing the fundamental operation of 
the instrument. 
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>V 1. INTRODUCTION 

SIM is designed as a space-based 10-m baseline Michelson optical interferometer operating in the visible waveband. 
This mission will open up many areas of astrophysics, via astrometry with unprecedented accuracy. Over a narrow 



field of view SIM is expected to achieve a mission accuracy of 1 /ias. In this mode SIM will search for planetary 
companions to nearby stars by detecting the astrometric "wobble" relative to a nearby (< 1°) reference star. In its 
wide-angle mode, SIM will be capable to provide a 4 /zas precision absolute position measurements of stars, with 
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parallaxes to comparable accuracy, at the end of a 5-year mission. The expected proper motion accuracy is around 
4 /^as/yr, corresponding to a transverse velocity of 10 m/s at a distance of 1 kpc. 1 

The SIM instrument does not directly measure the angular separation between stars, but the projection of each 
star direction vector onto the interferometer baseline by measuring the pathlength delay of starlight as it passes 
through the two arms of the interferometer. The delay measurement is made by a combination of internal metrology 
measurements to determine the distance the starlight travels through each arm, and a measurement of the central 
white light fringe to determine the point of equal pathlength. 

SIM surveys the sky in units called tiles. A tile is defined as a sequence of measured delays corresponding to 
multiple objects all made by a single baseline vector b and central pointing of the instrument - that is, all the mea- 
surements in a tile are from objects that are within a single astrometric FOR (field of regard of the instrument), which 
is approximately 15° x 15°. The existence of a single baseline vector insures that the system of equations developed 
from the observations to extract the astrometric parameters is not underdetermined. However, the collection of such 
a measurement set with a single interferometer is actually impossible, as the data collection on a sequence of objects 
takes finite time, over which both the baseline length and orientation do not remain constant. 

This paper describes the fundamental steps of how the on-board instrumentation of external metrology and 
auxiliary guide interferometers are used to reconstruct the baseline vector sufficiently accurately so that it can 
effectively be modeled as a single vector over the period of a tile observation. This process has been previously 
referred to as the regularization of the baseline. 2 The notion of the regularized baseline has been used extensively 
in a number of grid simulation studies that plan observation sequences, predict mission accuracy, and determine 
sensitivities to various instrument parameters. 2-4 

The process of reconstructing the baseline vector when implemented onboard in real-time is termed pathlength 
feedforward, and is a critical component to the operation of the interferometer. Because many of the astrometric 
targets will be very dim, it is not possible for the science interferometer to track the fringes and compensate for optical 
pathlength difference variations in real time using the dim target as the signal. As a result the fringes associated 
with the science target will be washed out due to uncontrolled motions of the instrument. The adopted solution in 
these cases is to use precise attitude information obtained from the two guide interferometers and construct a delay 
tracking signal that will be fed to the science interferometer's delay line in an open loop fashion. This aspect of the 
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Figure 1. SIM External metrology design 

baseline regularization process will be covered in some detail. An overview of the mission and several of the major 
subsystems of the instrument can be found in references 5-8 



SIM is designed to measure the pathlength delay between the two arms of the interferometer. The instantaneous 
delay value is given formally by the interferometer astrometric equation: 4 



where d is the external optical pathlength delay synthesized by a combination of internal metrology and white light 
fringe estimation, s is the normal to the wavefront of the starlight (the unit 3- vector to the observed object), b is 
the baseline 3-vector, A: is a so-called constant (or calibration) term that represents possible optical path differences 
between the light collected from the target object and the internal metrology, and r)(t) is the noise in the measurement. 

Because of limitations imposed by the optical throughput of the system and the brightness of the observed 
objects several seconds of integration time are necessary to bring the average white light fringe estimate to the 
required accuracy. Thus, the actual instrument measurement is the following: 



2. ASTROMETRY WITH SIM 



d(t) = (b(t) -s)+k + r](t) 



(1) 



(di) = • (6» +k + Vl 



(2) 



where (rfj) is the average measured external delay obtained from internal metrology measurements and white light 
fringe estimation, (b) is the average baseline vector over the period of the observation, and r]i is the measurement 
noise. 

The fundamental objective of the instrument is to make these delay measurements so that the astrometric 
parameters of position, proper motion, and parallax can be ascertained for the stars that are observed. Focusing just 
on the problem of estimating stellar positions, the astrometry problem is to determine the vector s in Eq. (2) above. 
However, all of the quantities above on the right must be treated as unknown because of the /xas level precision 
requirements of the instrument. (The Hipparcos catalogue has an accuracy on the order of several mas for stellar 
positions, while standard attitude determination and alignment systems on spacecraft can determine the baseline 
vector to the order of an arcsec.) SIM circumvents this difficulty by observing multiple stars, {s^}, within its field of 
regard so that the observation equations can be modeled as 

di = Si-b + k + rii, i = l,...,N. (3) 

The critical assumption here is that there is a single baseline vector, albeit unknown, that needs to be solved for 
as well as a single constant term. Thus, Eq. (3) consists of N equations with 2N + 4 unknowns. If a new baseline 
orientation is used to observe the same set of stars, an additional N equations are obtained, but only at the expense 
of 4 additional unknowns necessary to determine the new baseline vector and constant term. It is evident with more 
observations the system of equations eventually becomes overdetermined so that the stellar positions can be resolved. 
Applying this idea to a set of tiles that covers the entire celestial sphere is the kernel of SIM's strategy to perform 
wide angle astrometry. 

The observable in (2) is the average delay measurement made by the interferometer. Considerable analysis, 
simulation and technology development and validation has been devoted to the problems of precision white light 
fringe estimation and metrology that are used to synthesize this observable. 9-15 The focus of this paper is to show 
that the model equations Eq. (3) are actually consistent with the data that is collected by the instrument. To see 
where a potential conflict may arise, observe that because the baseline vector is time-varying, the average position 
over the integration period changes from observation to observation. The resulting model over a tile is 



(di) = Si ■ + k, i = l, 



N. 



(4) 



The difficulty now is that there is a different (unknown) average baseline for each star; violating the assumption 
made in (3). 

In principle SIM solves this problem by using guide interferometers and external metrology to track the baseline 
vector during the observations. The details of this process will be developed over the next two sections. Here will 
give the overview of how this is done, specifically with respect to overcoming the difficulty posed in (4). 

SIM uses two auxiliary guide interferometers that lock on bright "guide" stars (gi, §2), thereby keeping track of 
the directions to these stars, and hence also of the rigid-body motion of the instrument. The third interferometer 
switches between the science targets (si, s 2 , ■ ■ ■), measuring the projected angles between the targets and the interfer- 
ometer baseline vector. An external metrology system keeps track of the flexible-body motions of the instrument by 
measuring changes in the baseline vectors of the three interferometers in a local frame, or equivalently, determining 
their relative orientations. 

Because the three interferometer baselines are not collinear, to complete the characterization of the rigid body 
behavior of the instrument a third incrtial measurement is required. This measurement is termed the "roll" mea- 
surement. In the ideal case we show that the SIM instrumentation is sufficient in the sense that in the absence of 
measurement errors and a priori parameter errors, the collection of measurements made by the guide interferometers, 
the roll measurement, the external metrology measurements together with the a priori parameter data consisting of 
the positions of the guide stars, the initial guide and science baseline vectors in the local frame, uniquely determine 
the baseline vector of the science interferometer in inertial space. 

When this is the case the vector (bi) in (4) is known for all i. In reality there are both measurement errors and a 
priori parameter errors so it can never be assumed that the (bi) are known. For any parameter vector p, let b(t;p) be 
the estimate of the baseline vector using this parameter vector. In the absence of measurement error the sufficiency 
of the SIM instrumentation is embodied in the statement that b(t) = b(t;po) where po is the true parameter vector. 
Let p denote the a priori parameter vector and write 



Sb(t) = b{t;p )-b(t;p). 



(5) 



In Section 5 we show that for \p — po\ sufficiently small 

5b{t) = Sb° + e(t), (6) 

where 5bo is a constant vector (so long as the guide interferometers are locked on the guide stars) and \e(t)\ is a 
residual variation that contributes a delay error much less than the magnitude of the measurement noise r\i and can 
be ignored in the analysis. 

Now we may write the instantaneous delay equation as 

di(t) = Si-b(t;p) + Si-Sb{t) + k 

= • b(t; p) + s l - Sb° + ^ ■ e(t) + k. (7) 

Let Si° denote the a priori estimate of the position so that s*j = s*j° + Ssi, where 5si is the correction vector that is 
sought. Thus, after averaging 

(di(t)> - ■ (b(t;p)) = Ss* ■ (b(t;p)) + Si° ■ (5b°) + s*° • ( e (t)> + k. (8) 

The quantity on the left is termed the "regularized delay" and is the quantity that is synthesized from the SIM 
instrumentation. And since b(t;p) is known and the contribution of the term containing e(t) can be ignored, the 
unknowns are the corrections to the science star positions and the correction to the baseline vector which is now a 
single constant vector over the entire tile. In this sense the idealized model equations in (3) are correct. 
2.1. The ideal instrument measurements 

We will now get into the details of how (b) is obtained from the measurements and a priori data. For this purpose 
it suffices to treat SIM as a set of fiducials, Xi, ...,X N , e R 3 . (Here R 3 denotes Euclidean 3-space.) In the SIM 
reference design, shown in Figure 1, there are a total of six fiducials, and in a spacecraft local frame their coordinates 
are collected in the matrix X fld (with units in meters): 
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The i th row of X fid defines the coordinates of Xi. 



The science interferometer baseline vector b s is denned as 

b s = X lid (6,:)-X tid (l,:). (10) 

The auxiliary science baseline vector is defined by fiducials 3 and 4. The guide interferometer baseline vector b g is 
defined 

b g = X lid (5, :) - X fid (2, :), (11) 

and the "roll" vector is 

f=X fl d(2,:)-Xf id (l,:). (12) 

For any pair of vectors in R 3 , say Wi and Wj, the vector Wj will denote the difference Wij = W — Wj. Thus, 
for example b s — Xiq. Our interests center around the evolution of the fiducials Xi{t) over a time period to <t <T, 
where t = t denotes the beginning of an observation of a tile and t = T is the time of completion. The problem 
is solved using the on-board optical sensing systems that include the external metrology system, the guide star 
interferometers, and the roll estimator. The signals from these systems are briefly described next. 

In the SIM reference design relative distance measurements are made between each pair of fiducials except for the 
direct link connecting the active science interferometer fiducials. The observed variables associated with the external 
metrology system are 

t ij {t) = \X ij {t)\-\X ij (t )\, j > i,i = 1....5, excluding (i = 1 & j = 6). (13) 

These measurements are relative distance measurements, and (13) is valid for any choice of coordinate frame. Thus, 
b s is determined in the local frame from (10). The problem is to find b s in the inertial frame. This connection is 
made with the guide interferometers. 

SIM uses a pair of guide stars to produce two independent delay measurements per observation: 

d gk = (,9a • b gk ) + fc gA , (14) 

where <?a is the position vector to guide star A, with 4e (1,2). 

The roll estimator produces a "measurement" similar to the guide interferometers. We designate a guide telescope 
for use as part of the roll estimation scheme, and let n denote the line-of-sight vector of this telescope. Next we 



introduce a fiducial rigidly attached to the telescope which is also measured by the external metrology system. Let 
r denote the vector connecting this fiducial with the fiducial mounted on the chosen guide telescope. The rigidity 
assumption is that 

K = (n-T) (15) 

where k is constant for all values of n and r over the period of a tile observation. 
2.2. The logic of SIM astrometric observations 

Now we will describe how (13)-(15) are used to characterize the evolution of b s in the inertial frame. 

First let us see what can be learned from the observations in (13). Set X = (X\, Xn), and define F(X) as the 
function with components Fij(X) 

F ij {X) = \X ij {t)\-\X ij {t a )\, j>i; i = l, ...5, i + \ & j?6. (16) 

We seek the solution to the system of equations 

F{X)=i, e=(£ ij ). (17) 
The first thing to note is that if U is a rotation matrix acting on vectors in R 3 , and if £ G R 3 , then 

F(X) = F(Y), where X = (X u X N ), Y = (UXi + ...,UX N + |) (18) 

since 

IUX, + (- (UX 3 + 1)1 = \UX l3 \ = \X l3 \. (19) 

The final equality in (19) follows because U is orthogonal, and hence, preserves norms. Importantly the converse of 
(18) also holds: 

Property 1: If F(X) = F(Y) for some pair X and Y (in a small neighborhood N of X ), then there exists a 
rotation matrix U and 3-vector £ such that Y = UXi + £ for all i. 

This is the fundamental result which links the local and inertial frames. The linear justification of this principle 
is that since F(X) = F(Y) and to first order F(Y) = F(X) + F'(X)(Y - X), it follows that Y — X must be in the 
kernel of F'(X), which can be characterized as the rigid body motions of the fiducial system. It is not difficult to 
extend the linear argument to the full result. 



Let X denote the vector of fiducial positions in inertial coordinates. Then (in the absence of noise) X solves (17). 
If X loc is a solution to (17) computed in a local spacecraft coordinate frame, then Property 1 states that there is a 
rotation matrix U such that 

Xa = UX^ (20) 

for every pair of fiducials Xi and Xj. Thus the matrix U is the transformation between the local and inertial 
coordinate frames. And since the science baseline vector b s is known in local coordinates from external metrology 
measurements, the problem of determining 6 S in inertial coordinates is solved once we obtain U, viz. 

& s = C>fo s loc . (21) 

The equations for obtaining U are provided by the guide interferometer measurements and the roll estimator. 
The guide measurements may be written as 

rf g A = (ffA-?7& g loc ) + A :gA , (22) 

with 6g OC assumed known from external metrology data. The third equation needed to determine U is provided by 
the roll estimator, which has the form from (15) 

k= (n-f) = (n-Uf loc ), (23) 

where ft, k and r loc are all known. We now have 3 equations with which to determine U. (The set of orthogonal 
matrices live in a space of three dimensions, so three non-redundant equations are sufficient.) Note that the roll 
estimator equation has essentially the same form as the guide equation. The correspondences are that the guide 
measurement is replaced by a constant value using the rigidity assumption of the structure connecting the fiducials 
used for roll, the inertial position of the guide star is replaced with the inertial line of sight vector of the telescope, 
and the guide interferometer baseline vector is replaced with roll vector connecting the two fiducials. 

3. SOLVING THE BASELINE VECTOR EQUATIONS 

As described above, there are two components to the problem of determining b(t;p). The first part inverts the 1-D 
external metrology measurements into position vectors computed in the spacecraft local frame for all of the fiducials. 
The second part uses the guide interferometer measurements together with the roll estimator equation to determine 
the transformation between the local frame and the inertial frame. 



3.1. Inverting the external metrology measurements 

Because in general the system of external metrology equations Eq. (17) is overdetermined, the estimate of the fiducial 
positions Xi is derived from the solution to the nonlinear least squares problem 

min|F(X)-^| 2 , (24) 

where I is the vector with components £ij. In the case of the reference design, I is a 14-vector comprised of links 
between every pair of the six fiducials of the external metrology subsystem, except for the one pair that is a direct link 
between the fiducials of the active science interferometer baseline. Let F' denote the differential of the function F 
at X°. The rows of F' are constructed from the gradients of the functions Fij. These gradients are easily calculated 
analytically as 

(X? - X°) T - Xf) T 

VFy(Xo) = [0i, 3 (i-l), pfo_jfr)| ' °l,30-i-l), ~ |^0_jf0| ' °1.3(JV-j)], (25) 

where Stt is the zero matrix with s rows and t columns. Thus F' is a 14 x 18 matrix corresponding to the 14 
metrology measurements and 18 coordinates representing the 6 fiducial positions. Let F'^ denote the pseudoinverse 
of the differential of F at X°. (The pseudoinverse only needs to be calculated once at the beginning of the tile.) 
Then the iteration scheme beginning with X° = (X®, ...,X%), 

X k = X k ~ 1 +F'*[e-F(X k - 1 )], (26) 

can be shown to converge to the unique solution of (24) using a Gauss-Newton convergence argument. This solution 
is contained in a ball of radius 2\FU\ about X and an 0(\£\ 3 ) error is incurred if the iteration is stopped at k = 2; 
which is sufficient for both SIM requirements as the constraints on the magnitude of the flexible body motion of 
the fiducials leads to \£\ « 10 _5 m. The use of the pseudoinverse in (26) is equivalent to constraining certain linear 
combinations of fiducial positions to remove the rigid body motions of the system. 

3.2. The attitude equations 

The second part of the algorithm for estimating b{t) requires solving for the attitude. The three equations for 
obtaining the transformation, U, between the local and inertial frames are provided by the guide interferometer and 
the roll estimator equations. 



U can be parameterized in several ways. We will make use of the fact that there is a skew-symmetric matrix S 
such that U = cxp(5) so that U has the series expansion 



U = I + S + S 2 /2\ + S 3 /3l + .... 



(27) 



We will also use the 1-1 correspondence between the set of skew symmetric matrices and 3-vectors via the mapping 

LU — ► S(uj), 



S(<3) 



— 0>3 UJ2 

lo% — u>\ 

— L02 LOl 



, Q = (Wi, U> 2 , UJ 3 ). 



(28) 



Without loss of generality we may assume that U « / because of on-board attitude knowledge. The quadratic 
approximation to (22)-(23) using (27) and (28) is 



T ext co + Gext(w) = y, 



where, 





G(<3) 



and 



G(<3) = -2 



dgl - 51 • 



1g2 



K c x 5i 



(£ x 52, w x bl° c ) 



92 ■ bl 



(29) 

(30) 
(31) 

(32) 



It can be shown that the solution to (29) produces an error of 0(\uj\ 4 ) in determining U. (This is less than a 10~ 15 
rad error for 20 arcsec of motion of the instrument.) Without going into the details of the proof of this result here, 
we just remark that it hinges on the simple observation that for any skew symmetric matrix S, the distance from 
I + S + S 2 /2 to an orthogonal matrix is |0(5 4 )|. And this follows by noting that 

(J + S + S 2 /2)(J + S + S 2 /2) T = I + S + S T + SS T + S 2 /2+(S 2 f/2 + 

+ S(S 2 ) T /2 + S 2 S T /2 + S 2 (S 2 ) T /A = 
= I + S 4 /4, 



where / denotes the identity matrix and we have used S T 
orthogonal if XX T = I.) 



(33) 



S and (S 2 ) T = S 2 . (Recall that a matrix X is 



Rewriting (29) as 

Q = T-ly-T-\G^{i2), (34) 
the solution can be obtained by a fixed point iteration on the mapping defined above on the right: 

u) k = T-\y - T-\G^ t (fl*" 1 ), = T~\y- (35) 

with error estimate 

\Q k -u\<\T-^y\ k+ \ (36) 

A standard contraction mapping argument can be used to establish this result. The estimated baseline in inertial 
space is then realized as 

b s = bl° c + lu x ^ oc + ^x(ux ^ oc ), (37) 
where lr s oc is determined from (26), and uj is obtained from the iteration in (35). 

To simplify matters a little for the error analysis that will be performed in Section 5, we define ujj as the starting 
value in (35): 

<3/ = T~\y, (38) 

and bJii as the first iterate: 

^ // =^-T ext 1 G ex t(^). (39) 
Because our estimates will carry through second order, the science baseline estimate we use is 

b s = bl° c + Q u x ^ oc + \Qi x (<3/ x ^ oc ). (40) 

This second order approximation accommodates flexible fiducial motions on the order of tens of microns in the 
inversion of the external measurements and rigid body motions of 100 /^.rad in the attitude equations. These are 
both significantly larger than SIM requirements. As an example of these approximations, the plots in Figures 2 
and 3 show induced motion of this magnitude in the system. The x and z axis motion shown are due to large (20 
as) sinusoid motion of the baseline, while the y axis motion is due mostly to flexible body motion superimposed on 
the rigid body displacement. Figure 4 shows the error in the estimated baseline components using the quadratic 
model. It is seen that these errors are sub-picometer, as the error analysis predicts. These methods could of course 
accommodate a larger range of motion at the expense of utilizing more iterations in the solution of the external 



x 10" 




10 20 30 40 50 60 70 80 



90 100 



Figure 2. x and z baseline displacement vector components due to rigid body motion 

metrology and attitude equations. A mission level simulation capability that incorporates these equations and a 
number of instrument models is under development. 16 

4. A GEOMETRIC INTERPRETATION 

The solution to the set of non-linear equations that governs the evolution of the baseline vector has a geometric 
interpretation that will be discussed here. To facilitate this we introduce the baseline strain parameter e 



it) 



b(t) - b(t ) 

K*o) 



(41) 



which we expect to be ~ ^trad because the change in the baseline length will be on the order of 10 ^m. 

Again retaining terms through second order in u(t), we parameterize the instantaneous interferometer baseline, 
b(t), as follows: 



b(t) = b(t ) (l + e(t)) (n{t ) + [£3(t) x n(t )} + ^[£3(t) x p(t) x n(t )}] + 0(^ 3 )) , 



(42) 



where b(to), n(t ) are the initial baseline length and orientations of the guide interferometer; e(t) — Ab(t)/b(to) 
is the time-varying readings of the external metrology; uj(t) is the vector of small attitude changes in the baseline 
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Figure 3. y baseline displacement vector component due to flex body motion 
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Figure 4. Error in baseline vector estimate due to quadratic approximation 



orientation for the interferometer, uj(t ) = 0; and [a x b] denotes the cross vector product of vectors a and b. 
4.1. Science Interferometer Delay Equation 

The expression (42) allows us to write the delays for all three interferometers. Thus, to second order in uj one obtains 
the following expressions for the time-varying delays of the science and the two guide interferometers: 

d s (t) = k.(t) + b B fo)(l + e.(t))(3-{fkfa) + [W (43) 
d s (t) ee fcJi) + & g (i )(l + e g (i))(s A -{^ (44) 

Taking into account that, because of the flexible body motions, the rate of the attitude drifts of the guide 
interferometers is different from that of the science interferometer, we may write: 

tf g (t) = £J.(t) + [0 s (t) x C3.(t)] + 0(6l), (45) 

where 9 s (t) is the rate of temporal drift of the guide interferometer baseline's orientation relative to that of the science 
interferometer. Also, due to the flexible body motions, the interferometers may be misaligned at the beginning of 
the observations. We define the contribution of this misalignment vector, Wog(*o)i to the initial baseline orientation 
of the guide interferometer as below: 

n g (t ) = n s (t ) + Mo g (io) x n s {t )] + d(w„ 2 g )- (46) 

Subtracting the initial conditions from the time- varying delays tracked by the guide interferometers, Eq.(44), 
together with Eqs.(45),(46), define the instantaneous change in the science interferometer baseline orientation. How- 
ever, the obtained system of equations is underdetermined. This is why only two out of three components of the 
attitude drift vector cD s (t) may be determined this way. We call the undetermined component - the roll component 
- and denote it a s (t). For the determination of this component of the attitude vector SIM uses the roll estimator 
given by Eq.(23). 

In Section 3 we developed the iterations for solving the baseline estimation problem. Here we will present the 
final solution for the science interferometer in a form that is amenable to geometric interpretation. Let us first define 
the following notations: 

b(to) = b no — initial science baseline estimate (length, b — 6 B (to), and orientation, n = n s (to)); 



Sj, gi,ff2 — unit vectors for j th science star and that for the two guide stars; 
d s (t),dg(t) — instantaneous interferometric delays, of science and guide interfereometers; 

fc 0s — constant part of calibration term for the science interferometer at the beginning; 
Ak s (t), Afcg(t) — contribution of a temporal drifts in the calibration terms; 

t^og — initial misalignment of the guide interferometer baseline's orientation relative to that 
of the science interferometer; 
9g(t) — contribution of the temporal drift of the guide interferometer's orientation relative 

to that of the science interferometer (taken care by the external metrology); 
a s (t) — magnitude of roll of the science interferometer (determined from Eq.(23)). 

Given the solution of the attitude matrix, to s (t), and equation (43), and the presented notations, the instantaneous 
delay of the science interferometer is obtained in the following form: 

d s (t) = fcos + Afc s (i) + 6o(l + e s (i))|(no-Sj) + Az„(t)J, (47) 
with the feed-forward signal, Az ff (t), given by 

- \ (p^) [no x ft] P*{t) [no x ft] V-ff^j^ + 
z \ / [m- [91 x g 2 \) 

+ <*(*)((*>■[& >< ^) f:f X j}\ - {no-% x ^]) f:f X % 
V (n - [31X32J) {no-[9i x 32J) 



P g i(i) [n x ,g 2 ] -p g2 (i) [n x ft 



x 



/, Nr -, / -. \ r -,\ (no- [32 x 

(n -3i) [n xu>ogJ -(n o -0g(t)) [n x 31J ) — — —3 

V V 7 7 (no- [.9i x 52J) 

- ((n -ft) [n x w 0g ] - (n o -0g(i)) [n x 32]) fa ^ X Sj J 7 2 
v v 7 7 ("0- [.9i x ft]) 

+ o.(t) [no x (£ 0g - *,(*))] f [no x ft] ' [? X *}] > - [n x ft] fi^*^ ) + 

V (n - [31X32J) (n - [31 x 32J)/ 

+ o(jV(t);co 2 ;el(t):9l(t)), (48) 



where the quantity p s (t) characterizes the pathlength feed- forward signal (instrumental drifts) and has the form: 



A variation in the average calibration term, k, within a tile has implications on the overall astrometric accu- 
racy. Such inner-tile k variation translates to a distortion of the relative geometry for the objects in the tile. The 
implications of inner-tile k variation on the grid reduction accuracy is current topic of study. 
4.2. Geometric Interpretation 

The expressions for the science delay (47)-(49) obtained above have a clear geometric interpretation. A discussion 
of the significance of the different terms contributing to the science delay is taken up below. 

The first two terms in (47) are the constant delay offset and the total temporal variation in this delay since 
the beginning of the tile observation. In the current SIM reference design, 5 the challenge is to keep the temporal 
variation in this term small, such that its total contribution to the science delay will be negligible (on the order 
of tens of picometers). The terms in the big parentheses are scaled with the change of the baseline's length, e(t); 
the magnitude of this scaling function is on the order of 1 ^rad, (c.f. (41)). The first term in the parenthesis 
is the usual form for the interferometric delay (b ■ s) , which represents the constant delay at the beginning of the 
observations of the tile. The second term in this parenthesis is the feed- forward signal, Az ff (t), which is given in (48). 
The first two terms in Eq. (48) are the first- and second-order feed-forward signals that update the instantaneous 
attitude of the science interferometer's baseline vector as measured in the local inertial frame. The third term in this 
expression is the second order contribution introduced by the temporal changes in the baseline lengths of the two 
guide interferometers (e g , e g ) which are coupled to their initial misalignments (w g, Wo g ) with respect to the science 
interferometer's baseline vector. The next term in Eq. (48) is due to the interaction between the first order feed- 
forward signal and the rigid-body motions of the entire extended structure. It's instantaneous contribution is given 
by the accumulated temporal drifts in the guide interferometers baselines orientation over the time (0 g , g ). The last 
term - is the second order contribution of the science interferometer's baseline roll vector (given by a s ) to the total 
science interferometer's delay. 

One of the results that immediately follows from the presented analysis is the criteria for choosing the positions 
of the guide stars in the tile. Indeed, one may see that all the terms responsible for the feed-forward signal have 
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Figure 5. Error in baseline vector estimate due to guide star position error 

the quantity (n - [gi x g 2 ]) in the denominator, see Eq. (48). This quantity is proportional to the volume of the 
tetrahedron constructed from the three vectors, namely the unit vector of baseline orientation ri and the directions 
to the two guide stars <?i and <f 2 - As this quantity appears in the denominator, it is preferable to have it as large 
as possible, thus minimizing the first and second order terms. In addition, by minimizing their contributions, one 
minimizes the contribution of guide star position errors to the estimate of science delay. The SIM reference design 
is near optimal in the sense that the guide stars are placed 90° apart, with each guide star nearly orthogonal to the 
interferometer baseline vectors. 

5. BASELINE REGULARIZATION CONCEPT 

When all the parameters used to determine the baseline vector via the prescription developed in Section 3 are without 
error, then the estimate of the baseline vector is constructed without error. But this can never be the case, and in 
Section 2 we showed that it is sufficient that the error be constant. This is the linchpin idea for the operation of the 
instrument. 

In Figure 5 the effect of 20 mas errors in the guide star directions are introduced in determining the baseline 
vector. The guide star position error results in an error in the estimate of the attitude motion of the system, leading 
in turn to an error in the baseline vector estimate. The x axis error is observed to be on the order of .0.77 ^m. 



However, observe also that this value is highly stable with a peak to peak range of approximately 10 pm. Thus 
although b s (t) — b s {t;p) exhibits micron class error, the important aspect is that the error is very stable as required 
in (6). 

This example provides corroboration of the notion of the regularization concept introduced in Section 2. Now we 
proceed to give a more rigorous mathematical justification of regularization in the face of a priori parameter errors, 
such as the guide star position error discussed above. 

Returning to the prescription for determining b(t; p) in Section 3, a priori values must be provided for the following 
parameters: 

Guide Interferometer Parameters: 

gi = guide star position vectors (i = 1, 2) 

fcj = guide interferometer constant terms (i = 1, 2) 

Roll Estimator Parameters: 

k = inner-product of roll vector r with guide telescope unit line-of-sight vector 
Optical Truss Parameters: 

X a = initial vector of all fiducial positions in a s/c local frame 
The totality of these parameters constitute the parameter vector p, 

P = (<7i,<?2, c 1 ,c 2 ,k,X°) (50) 
that is used to define the baseline estimate. In what follows we will assume the true parameter vector is 

Po = + 5gi,g 2 + 5g 2 ,c 1 ,c 2 , k + Sk,X° + SX°), (51) 
while the nominal parameter vector is 

p=(g u g 2 ,0,0,K,X°). (52) 



Hence the error parameter vector is the difference 



Po-p = (Sgi,Sg 2 , ci, c 2 , 5k, 6X°). (53) 

We will now prove a "weak" version of regularization that asserts to first order (the meaning of which is made 
more precise below) the difference between the baseline estimate and the true baseline vector is constant. Some 
comments on the second order errors that are ignored in this analysis are made at the end of this section. 

The following nomenclature will be used in the course of establishing the regularization result. 

Zeroth order: A quantity expressed in meters or radians (e.g. the length of an interferometer baseline vector.) 

First order: A quantity expressed in terms of um or /irad, or a product of a first order term with a zeroth order 
term (e.g., the magnitude of the elastic deformation of the optical truss, s/c attitude stability, guide star position 
error.) 

Second order: A product of first order quantities, or a product of a second order term with a zeroth order term. 

Let 6b s denote the initial error in the science baseline vector stemming from the fiducial error term 5X°. Also let 
a)/, (on denote the attitude rotation vectors in (38)-(39) obtained using the parameter vector p instead of the true 
vector p . Then by taking differences we obtain 

b s (t;p )-b s (t;p) = Sb° s + cJ/ x Sb° s + {w n - wjj) x ^ oc + i(cJ/ - d)j) x (<3j x ^ oc ) + 

+ ^ujj x ((<& - Cjj) x ^° c ) + - uj!) x [(tfj - wj) x ^ oc ]. (54) 

Retaining first order terms according to the definitions above, 

b.(t,po) - b s (t;p) = 6b° + (aJjj - u n ) x ^ oc . (55) 

This can be further simplified to 

b,(t,po) - b s (t;p) = Sb? + (tfj - ui) x ^ oc (56) 

because to first order Cjj = tin and wj = Uu since the iteration in (39) is a second order correction (see definition 
of G in (31)). Once we show that (w/ — u)j) is constant to first order, we will be done since 6^ oc is constant to zeroth 
order (its variations are due to elastic motions of the instrument). 



With these definitions the true (second order) attitude equation is given by 



T ext w + Gext(^) = y + Sy, 



(57) 



while the attitude equation using the parameter vector p is 



where y is defined in (32) and 



with 



T ext uj + G ext (w) = y, 
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Examining these terms it is seen that to first order 5y is a constant. 

The equations (57)-(58) are solved iteratively to obtain the true and estimated values of uj. Thus we have 



(62) 



and 



w/ = T-\y. 



(63) 



Hence, 



ui - (bi = T ex lSy. 



(64) 



Because T ext is a constant matrix to first order and Sy is constant to first order, it follows that ljj — Qji is also constant 
to first order; as was to be shown. 

It is possible to develop very specific forms for the second order error terms due to these initialization errors, but 
this is beyond the scope of the present paper. We remark that these errors do drive certain requirements in the SIM 
astrometric error budget. For example there is a relatively strong coupling between the initial fiducial error 5X° and 
attitude motion that produces a second order delay error of the form 



d(t) = [u Jl (t)x(st-Sb g l S^)]-s. 

H oc \ 



(65) 




Figure 6. Delay error due to absolute metrology and ACS errors. 

In Figure 6 we plot this delay regularization error for a "typical" case arising from a single gauge absolute metrology 
error of 5 ^m in determining the fiducial positions and a 2 arcsec single axis rigid body motion generated from a 
random rotation direction with a 10 /ir&d amplitude sinusoidal motion. The plot contains the associated delay error 
for a typical star within the field of regard of the science interferometer. If the abscissa is the time-axis in seconds, 
and say a 15 sec observation of the star is made, then the delay error of the observation would be the average value 
over that particular 15 sec period. The observation period could be at the extremes of the motion, producing a delay 
error of ±150 pm, or it could be centered around zero producing zero delay error. In fact any delay error could be 
attained within this envelope. Requirements are set to keep this envelope small. 

6. SUMMARY, CONCLUSIONS AND FUTURE PLANS 

SIM astrophysical science is extracted from a model set of equations that relate the measured optical pathlength delay 
to the projection of the interferometer baseline vector onto the star direction vector. These equations presuppose 
that the interferometer baseline is fixed in inertial space; which it is not. The main objective of this paper has been 
to introduce the reader to the concepts and the instrumental logic of the SIM astrometric observations, especially 
as they relate to the fundamental operation of baseline regularization that "fixes" the interferometer in inertial 
space. Mathematical arguments were presented to establish this fundamental principle and a precise definition of 



rcgularization was given. The underlying nonlinear system of equations that is the basis for regularization was 
derived and numerical methods to solve them were obtained. A simulation was also developed incorporating the 
numerical processing methods of the instrument observables and the results were shown to conform with the theory. 

Beyond demonstrating the SIM proof of principle, the regularization equations presented here also form the 
kernel of the extensive instrument error budget. The linearized version of these equations are used to determine the 
propagation of noise from external metrology measurements and guide star delay measurements to the science delay 
error. This error is determined by the geometry of the optical truss together with its orientation in inertial space 
with respect to the guide star and target star positions. The propagation factors are used to set requirements on 
the integration time of observations, single gauge metrology error, etc. The regularization equations also reveal the 
existence of a number of second order errors that arise in the form of products of fiducial motion (elastic and rigid 
body) and initial parameter error. Examples of errors of this type were given. 

Current work focuses on the mechanisms and effects of variation of the "constant" terms in the astrometric delay 
equations. Nominally the appearance of the constant term compensates for the lack of a precise internal metrology 
gauge that measures the absolute distance between the interferometer aperture fiducials to the beam combiner. 
However, there are a number of instrument errors that are collected into this single term. For example as the 
interferometer observes stars within its field of regard several optical elements must be translated and rotated. These 
induce non-trivial diffraction effects, metrology gauge error due to imperfect corner cubes, reflection phase errors 
because of a changing angle of incidence of the interrogating metrology beams, and others. Each of these effects 
must be played through the delay regularization equations to ascertain their ultimate effect on delay error. 
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